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Abstract 

We study the behavior of quark and diquark condensates in dense quark matter 
under the influence of a gravitational field adopting as a simple model the static 
D— dimensional Einstein Universe. Calculations are performed in the framework of 
the extended Nambu-Jona-Lasinio model at finite temperature and quark density on 
the basis of the thermodynamic potential and the gap equations. Quark and diquark 
condensates as functions of the chemical potential and temperature at different 
values of the curvature have been studied. Phase portraits of the system have been 
constructed. 



1 Introduction 

As is well known, non-perturbative effects in low energy QCD can be studied with the 
use of various effective models with four-fermion interaction, such as the instanton model 
PP or the Nambu-Jona-Lasinio model (NJL) [2]. 

Effective field theories with four-fermion interaction of the Nambu - Jona-Lasinio 
type, which incorporate the phenomenon of dynamical chiral symmetry breaking, are 
quite useful in describing the physics of light mesons (see e.g. [31 HI [5] and references 
therein) and diquarks [SI El- Moreover, on their basis one can study the effects of various 
external conditions, like temperature and chemical potential [H] , and consider the influence 
of external electromagnetic and gravitational fields. In particular, the role of the NJL 
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approach increases, when detailed numerical lattice calculations are not yet admissible 
in QCD with nonzero baryon density and in the presence of external gauge fields [HI [TD1 
[TTl [T2l [T3] . Moreover, the N JL model has also applications in the nuclear physical and 
astrophysical researches (e.g., neutron stars) as concerns the quark matter [14]. 

It was proposed more than twenty years ago [H2, [161 Ej that at high baryon densities 
a colored diquark condensate < qq > might appear. In analogy with ordinary supercon- 
ductivity, this effect was called color superconductivity (CSC). In particular, the CSC 
phenomenon was investigated in the framework of the one-gluon exchange approxima- 
tion in QCD [18], where the formation of diquarks (colored Cooper pairs of quarks) is 
predicted self-consistently at extremely high values of the chemical potential fi > 10 8 
MeV [19] . One should also mention the use of the Bethe-Salpeter equation to obtain the 
diquark masses in the CSC phase of cold dense QCD at asymptotically large values of the 
chemical potential [2D]- Unfortunately, the corresponding high baryon densities are not 
observable in nature and not accessible in experiments. In order to study the problem 
at lower values of /z, various effective theories for low energy QCD, such as the instanton 
model p] and the Nambu— Jona-Lasinio (NJL) model|2] can be employed. The possibil- 
ity for the existence of the CSC phase in the region of moderate densities was recently 
proved (see, e.g., the papers [U [2H [221 [23]) as well as the review articles [2H 122] and 
references therein). In these papers it was shown that the diquark condensate < qq > 
can appear already at a rather moderate baryon density (/1 ~ 400 MeV). The conditions 
favorable for this condensate to be formed can possibly exist in the cores of cold neu- 
tron stars. Since quark Cooper pairing occurs in the color anti-triplet channel, a nonzero 
value of < qq > means that, apart from the electromagnetic U(l) symmetry, the color 
SU C (3) symmetry should be spontaneously broken inside the CSC phase as well. In the 
framework of NJL models the CSC phase formation has generally been considered as a 
dynamical competition between diquark < qq > and usual quark- ant iquark condensation 
< qq >. Special attention has been paid to the catalyzing influence of the vacuum color 
fields on the condensation of diquarks [T3] (see also, [2D])- Moreover, mesons and diquarks 
in the color neutral superconducting phase of dense cold quark matter were studied in 
[27] and Nambu-Goldstone bosons in the color-asymmetric superconducting phase were 
considered in the framework of the Nambu-Jona-Lasinio-type model in [28J. 

In several papers, in the framework of the NJL model, the influence of a gravitational 
field on the dynamical chiral symmetry breaking (D%SB) due to the creation of a finite 
quark condensate < qq > has been investigated at zero values of temperature and chemical 
potential [291 EDI EH E2]. The study of the combined influences of curvature and tem- 
perature has been performed in [33]. Recently, the dynamical chiral symmetry breaking 
and its restoration for a uniformly accelerated observer due to the thermalization effect of 
acceleration was studied in [34J at zero chemical potential. Further investigations of the 
influence of the Unruh temperature on the phase transitions in dense quark matter with a 
finite chemical potential, and especially on the restoration of the broken color symmetry in 
CSC were made in [35]. One of the widely used methods of accounting for the gravitation 
is based on the expansion of the fermion propagator in powers of small curvature [36l [37] . 
For instance, in [3H], the three-dimensional Gross- Neveu model in a space-time with a 
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weakly curved two-dimensional surface was investigated, using an effective potential at 
finite curvature and nonzero chemical potential. In paper [2H], this approximation was 
used in considering the D^SB at non-vanishing temperature and chemical potential. It 
should, however, be mentioned that near the phase transition point, one cannot consider 
the critical curvature R c to be small and therefore exact solutions with finite values of R c 
should be used. This kind of solution with consideration for the chemical potential and 
temperature in the gravitational background of a static Einstein universe has recently 
been considered in [40J. In all these papers it was demonstrated that chiral symmetry 
is restored at high values of the space curvature. There arises the interesting question, 
whether similar effects can occur in the study of CSC under the influence of a gravitational 
field. 

In the present paper, we study the influence of temperature, chemical potential and a 
gravitational field on the creation of the CSC state in quark matter. As is well known, one 
of the possible models for the expanding universe is the closed Friedman model. However, 
since the formation of condensates is expected to take place considerably faster than 
the rate of the Universe expansion, its radius in our calculations can be considered as 
constant. Therefore, as a simple model of gravitation with a given background space-time 
we assume a static Einstein universe. In this case, we will derive an analytical expression 
for the effective potential of the model by using the mean field approximation, with all 
the information concerning the chiral and color condensates. As is known, in a flat space- 
time, the chemical potential and temperature lead to fairly different effects, i.e., at large 
chemical potential, a first order phase transition takes place, while high temperatures 
lead to the second order phase transition. The aim of this paper is to study the phase 
structure of the quark matter, in particular CSC, in a highly curved space-time in the 
model of the static Einstein universe. 

2 The extended NJL model in curved space-time 

In D-dimensional curved space-time with signature (+, — , ■), the line element is 

written as 



The gamma- matrices 7^, metric g^ v and the vielbein e^, as well as the definitions of 
the spinor covariant derivative V„ and spin connection uj^ b are given by the following 
relations [SZHH]: 



ds 2 = •q & ie & e b v dx tM dx v . 




diag(l, -1,-1, -!,•••), 



(1) 



(2) 
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Here, the index a refers to the flat tangent space defined by the vielbein at space-time 
point x, and the 7 a (a = 0, 1, 2, 3, • • •) are the usual Dirac gamma-matrices of Minkowski 
space-time. Moreover 75, is defined, as usual (see, e.g., [37J, H21 HH] ) , i.e. to be the same 
as in flat space-time and thus independent of space-time variables. In the following we 
shall consider the static Einstein universe with the line element 

ds 2 = dt 2 — gij{x)dx % dx^ (i, j = 1, . . . , D — 1). (3) 

A conventional model that demonstrates the formation of the color superconducting 
phase is the extended NJL model with up- and down-quarks. This model may be consid- 
ered as the low energy limit of QCD. For the color group SU(3) C its Lagrangian takes the 
form0 



C = q [i^Vft + q+St (qq) + (qrfrq) 



(4) 



+ ^ [iq c ee b "/ 5 q] [iqee b -f 5 q c ] . 

Here, /i is the quark chemical potential, q c = Ctf ,q c = q l C are charge-conjugated 
bispinors (t stands for the transposition operation). The charge conjugation operation 
is defined, as usual (see, e.g., f37j), with the help of the operator C — ry 2 7°, where the 
flat-space matrices 7 2 and 7 are used. 

The quark field q = q lol is a doublet of flavors and triplet of colors with indices 
i = 1,2; a = 1, 2, 3. Moreover, f = (r 1 , r 2 , r 3 ) denote Pauli matrices in the flavor space; 
(e) lk = e lk 1 (e b ) a P = e af3b are the totally antisymmetric tensors in the flavor and color 
spaces, respectively. 

Next, by applying the usual bosonization procedure, we obtain the linearized version 
of the Lagrangian (jlj) with collective boson fields a, 7? and A, 

£ =q[i^V^ + ^]q-q(o + ir > TTi)q-^- i (a 2 + T? 2 )- 

-^-A* b A b - A* b [iq* Cee b ^q] - A b [iqee b ^ Cq l ] . 1 } 

The Lagrangians (jlj) and ([5]) are equivalent, as can be seen by using the equations of 

motion for the boson fields 

\ h G2 . t „ k 5 G\ _ _ G\ _, h _ 

A = — —iqCse-iq a = — —qq, n = — — 527 rq. 

The fields a and 7? are color singlets, and A b is a color anti-triplet and flavor singlet. 
Therefore, if < a > 7^ 0, the chiral symmetry is broken dynamically, and if < A 5 > 7^ 0, 
the color symmetry is broken. 

3 Effective action 

Following [13], we obtain the partition function of the theory 

Z = J [dq] [dq] [da] [dn] [dA* b ] [dA b ] exp j i J d D x^~g I } , (6) 

1 In this work, we shall consider the particular realistic case of N c ~ 3. The general case of arbitrary 
N c is discussed in [1] with relation to the effective hadron theory of QCD. In what follows, however, we 
shall keep the general notation N c for convenience. 
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where g = det | g^. In what follows, it is convenient to consider the effective action for 
boson fields, which is expressed through the integral over quark fields 



exp { iS e{i (a, vr, A b , A* b ) } = N' J [dq] [dq] expj i J d D x^~gt 



(7) 



where 



* b \ 



d D x, 



-g 



3(a 2 + 7f 2 ) 



+ 



3A b A* 



+ S q , 



2G\ G2 

N' is a normalization constant. The quark contribution to the partition function is given 
by the expression 



Z q = e lb " 



N' J[dq][dq]esxp[ i J d 



id 



q (ry^V^ — a — ij 5 nr + ^7°) q — 



-q (iA b ee b j 5 C) q 4 - q l (iA* b Cee b ^) q 



(9) 



In the mean field approximation, the fields a, tt, A b , A* b can be replaced by their ground- 
state averages: <cr>, < n >, < A b > and < A* 6 >, respectively. Let us choose the 
following ground state of our model: <A 1 > = <A 2 > = <7r> = and < a >, 
< A 3 > 7^ 0, denoted by letters a, A. Evidently, this choice breaks the color symmetry 
down to the residual group SU C (2). 

Let us find the effective potential of the model with the global minimum point that 
will determine the quantities a and A. By definition S e g = — V e g J d D x^—g, where 



3a 2 3 A A* ~ 
V eS = — + —=— + Vdf ; 



2d 



Go 



v 



d D Xy / —g. 



(10) 



In virtue of the assumed vacuum structure (< A 1,2 >= 0), the functional integral for S q 
in © factorizes into two parts 



expiS^cr, A) = N' J[dq 3 ][dq 3 ] exp (i J dPxyf^g q 3 Dq 3 ) x 



x / [dQ] [dQ] exp (if d D x^g QDQ + QMQ 1 + Q l MQ 



where q 3 is the quark field of color 3, and Q = (qi, q^) 1 is the doublet of quarks of colors 
1,2. The following notations are also introduced 



D = i^Vf, -a + w c 



M 



y-,5 



-iA* Cee-y 



M 



-iAeej 5 C. 



(12) 



Here e = (e) 0/3 denotes an antisymmetric tensor in the color subspace corresponding to 
the SU C {2) group. 

Clearly, integration over q 3 in ( ITTi) leads to DetD. Let us next define \I/ = (Q*, Q) 
and introduce a matrix operator 

2M, -D* 
D, 2M 



Then the Gaussian integral in fill I) can be written in a compact form 

J [^]e5/ dD ^^* = Det i/2 Z . 



(13) 



Using the general formula 
A B 



Det 



x B A 

we obtain 

exp (iS q (a,A)) 



\S A ) = Det ^~ BB + BAB ~ lA ] = Det i AA ~ ABA ~ l B] 



N' Det ID) Det 1 / 2 



+ MD t M~ l D 



TV' Det 
x Det 1 /* 



(W — a + /i7 C 



x 



(14) 



4|A| 2 + (-iV - tr + //7°)(iV - a + /i7 c 



Here, the first determinant is over spinor, flavor and coordinate spaces, and the second 
one is over the two-dimensional color space as well, and V = 7 M V M . 

4 Evaluation of determinants in the static Einstein 
universe 

Let us first calculate the contribution of the determinant 



In Det 



(iV -a + /ry°) 



\ In Det 
I In Det 



(iV - cr + /i7°) 



+ i In Det 



7 5 (iV - a + yU7°)7 5 



In our case, r is evidently equal to zero, and g^V^Vj, = — V 2 , where V 2 = g^VjV 
is the spinor Laplacian. Then 



2 In Det 



(iV - a + /i7 C 



In Det 



In ( 1X51) we used the identity 



i 



(15) 



(16) 



where R is the scalar curvature of the space-time. Next, consider the contribution of the 
second determinant 



Det 



4|A| 2 + (-zV - a + //7°)(zV - a + /i7 c 



Det [4|A| 2 + 7 M 7"V M V, + a 2 + /i 2 - 2/ia 7 + 2i/ryVV*] 



(17) 
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In analogy with the case of a flat Minkowski space-time, consider the Hamiltonian of a 
massless particle 

H = aft (18) 

where a k = 7°7 fc , and (p) k = — zV^, k = 1, . . . , D — 1. One can easily demonstrate that 
in the case of the static metric (ED 



H 2 



-V 2 + -R. 



(19) 



With consideration of ffTB"]) and (fT9"j) . formula (fTTI) can be rewritten in the form 



Dot 



4|A| 2 + (-iV - a + Ai7°)(«V - cr + /X7 C 



Det 



4IAI 



<r 2 + fi 2 -p 2 + H 2 -2fi (# + a 7 

where po = ^o- Finally, let us introduce the Hamiltonian of a massive particle 

H = ap + <J7°, 

where p = — iV. It is related to the massless operator as follows: 

fe = h 2 

Then, (1171) is rewritten in the form 



a 



(20) 



(21) 



:22i 



Det 



4|A| 2 + (-iV - a + //7°)(«V - a + /i7° 



Det 



4|A| 2 + /i 2 -p 2 + ft 2 -2/iH 



(23) 



Thus, in order to calculate the determinants defining the effective action, one has to 
find the eigenvalues of the operators H and 7i. In particular, they may be found exactly 
for the case of the static D-dimensional Einstein universe. They are expressed through 
the corresponding eigenvalues of the Dirac operator on the sphere S 15-1 [4"3l 14*4"] . The line 
element 

ds 2 = dt 2 - a 2 (d6 2 + sin 2 6d^l D _ 2 ) (24) 

gives the global topology K. (g> S D_1 of the universe, where a is the radius of the universe, 
related to the scalar curvature by the relation R = (D — 1)(D — 2)a~ 2 . The volume of 
the universe is determined by the formula 



V(a) 



2n D / 2 a D - 1 

r(f) 



(25) 



Let us denote the absolute values of the eigenvalues of the operators ffT8l) and fl2TT) by 
lo\ and Ei respectively. Then, according to [131 HI] one has 



±ui4> h uji 

±E l i/> l , Ei 



/ = 0,1,2... . 



(26) 



The degeneracies of uj\ and E\ are equal to 



2lWlT{D + I - 1) 

dl ~ nv(D-i) ' (27) 



where [x] is the integer part of x. 



5 Thermodynamic potential and gap equations 

Next, let us take the complete system of eigenf unctions ^(t, x) = e~ ipot ipi( basis. 
Then the operators in the two determinants are diagonal and the latter are reduced to 
products of corresponding eigenvalues Ef — (po — fi) 2 and 4| A| 2 — p§ + (Ei ± fi) 2 . Evidently, 
there appears a gap in the spectrum of the second operator proportional to |A|, and as 
a result, the diquark condensate is formed. Using the relation detO = exp(trlnO) and 
following the standard procedure, one obtains the following expression for the quark 
contribution to the effective potential 

~ OO f 

V eS = ~% =i^ftrT l Ed l {H(E l ±ri 2 -p 2 }+ f 

1=0 ± L (28) 

+21n[4|A| 2 -p 2 + (^± / i) 2 ]}, 



where v = Vt and V is the volume of the universe, defined in ( f25j) . t is the time interval, 
and Nf denotes the number of flavors Nf = 2. 

In the case of finite temperature T = 1/(3 > 0, the following substitutions in ( |28i) 
should be made: 



where u n is the Matsubara frequency. Then the quark contribution to the effective po- 
tential (|28|) becomes the thermodynamic potential Q g 

n=oo oo f 

n, = E EE^{in[(^±^) 2 + ^] + 

P n=-ooZ=0 ± L (29) 

+21n[4|A| 2 + cu 2 + (^±/i) 2 ]}. 

The summation over the Matsubara frequencies can now be performed with the use of 
the formula 



p 2 

^lnK + p 2 )=^/ rf x 2 -^ + ^ln(. 2 + l) 



(30) 
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where, according to (I2"9"j) . p denotes a/ {Ei ± p) 2 or y4j A| 2 + (2^ ± /i) 2 , respectively. In 
the following, we shall neglect the last term in (I30p . since it is independent of a and A, 
and the first term can be rewritten as 

p 2 

Yl / d%2 ~2 2 = 21ncosh(p/3/2) + const = p(3 + 2\n (l + e~ p/3 ) + const. (31) 

J W n + 3; 

1//3 2 

Then the final form of the thermodynamic potential looks like 



(32) 



fi(<x,A) = 3(^ + ^ 

oo 

— ^r(N c -2)Y,di {Ei + Tin (l + e -W-")) + T In (l + e -^+A*)) } 

CO 

E di{ - ^) 2 + 4 I A I 2 + + Z^) 2 + 4 I A I 2 + 

z=o L 

+2Tln (l + e^V(^W) + 2Tln (i +e ^V(WW) }. 

In order to find the correspondence of the result obtained with the case of the flat 
space-time 0, we consider the zero curvature limit, i.e., R — > (a — > oo), of the expression 
for the thermodynamic potential. In this limit 

i? — ► 0, a — ► oo, I — ► oo, cj/ — > /a -1 , — > J dl = a J dui = a J dp 

(where we have changed the variables uj\ — > p) and, since 

I 1 



r(z + i) 

the degeneracy is equal to 

2KWV2]£p + / _ i) ^ 2 t( D+1 )/ 2 ] D _ 2 

' ~ r(i + i)rp-i) ~ (d-2)\ 

Then, taking expression (1251) for the volume V(a) into consideration, we finally obtain 

-2N f j -f^ {E p + Tin (1 + e -^-/*)) + Tin (l + e ~^ E ^) } - 

~ 2N f I (0{ V TE p - ^) 2 + 4 | A| 2 + y/(E p + p) 2 + 4| AF+ (33) 
+2Tln (l + e -/V(^) 2 +4|Apj + 2Tln ^ + e -^{E p+ rf+^ j 



2 Due to different topologies of the flat space-time and the Einstein universe, the expressions for their 
line elements do not relate to each other by any simple limiting procedure. 
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This equation coincides with the corresponding result for the flat space-time (see, e.g., 

[211 1221123). 

Now, imposing the condition on the effective potential, f2(0, 0) = 0, we should subtract 
a corresponding constant from it. The thermodynamic potential, normalized in this way, 
is still divergent at large I, and hence, we should introduce a (soft) cutoff in the summation 
over I by means of the multiplier e~ Wi ' , where A is the cutoff parameter. 

As is known, in the flat space-time, the regularization cutoff constant A can be deter- 
mined from the experimental results in measurements of the pion mass or the pion decay 
constant. However, in the case of a curved space-time, there are no experiments of this 
kind, and in principle we cannot determine the value of the cutoff parameter. Therefore, 
similar to [ID], we shall concentrate primarily on a qualitative discussion of the results 
obtained, concerning the phase structure of the universe. 

To this end, we shall divide all the dimensional quantities that enter the thermo- 
dynamic potential by the corresponding power of A to make them dimensionless, i.e., 
n/A D , a/A, A/A, A D - 2 G lt2 , A^V, R/A 2 , T/A, /i/A, u t /A, and denote them us- 
ing the same letters as before: Q, a, A, G 1)2 , V, R, T, fi, ui. Then the regularized 
thermodynamic potential is written as 

^(a,A)=3(^ + ^)- 

OO 

-y-{N c ~ 2) E e ~ Uld i { E i + Tin (1 + e-^*-")) + Tin (l + e -^ El+ ^) } - 

(34) 



N 



1=0 

oo 



/ E e-^dA v/(^-/i) 2 + 4|A| 2 + v /(^ + / u) 2 + 4|A| 2 + 
i=o 1 

+2T In (l + e-/V(^) 2 +4|Apj + 2T In (l + e -Py/W+tf+m^ } . 

In the following Section, we shall perform the numerical calculation of the points of the 
global minimum of the finite regularized thermodynamic potential f2 reg (<7, A) — f2 reg (0, 0) 
(they should of course coincide with the minima of the potential f2 reg (cr, A)), and with 
the use of them, consider phase transitions in the Einstein universe. To this end, we have 
to consider the following gap equations: 

<9Q reg <9fi reg 



6 Phase transitions 

As it is seen from ( 132*1) . the thermodynamic potential depends on two independent coupling 
constants G\ and G2. In what follows, we shall fix the constant G2, similarly to what has 
been done in the flat case [221 "13] , by using the relation 

G 2 = \g x . (36) 

o 

For numerical estimates, let us further choose the constant G\ such that the chiral and/or 
color symmetries were completely broken. For comparison, we shall consider two cases of 
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the choice of the constant: 1) G% = 10 (this value formally corresponds to that of |40j ) 
and 2) d = 20 @. Moreover, let us now limit ourselves to the investigation of the case 
D = 4 only. 

Let us first consider phase transitions in the case of zero temperature. Then, the 
thermodynamic potential takes the form 

n«*(a,A) = 3(^ + ^)- 

oo 

-%-(N c - 2) £ e-** {E t + (/i - EfiOQA -E t )}- (37) 
i=o v ' 

E e- u 'd, { - )") 2 + 4 I A I 2 + + V) 2 + 4|A|4 . 

2=0 L J 




Figure 1: Condensates a and A as functions of \x for i? = 3, Gi = 10, T = (all quantities 
are given in units of A). 



Figs. [H [21 E] show the behavior of the point of the global minimum (cr, A), determined 
from the gap equations (IHol) for the potential (13TI) . as a function of /i at different values 
of R = 3, R = 12 and T = 0, for d = 10 and d = 20. 

Using the formulas for the thermodynamic potential f2 reg , we can calculate the fermion 
number density in the system 

<9fp g 

n = — at fixed T. (38) 

The function n(/i) at fixed R is depicted in Fig. HI The step-wise character of the curve is 
explained by the discreteness of the energy levels and by the presence of the step-function 
in (13 7p . It is seen from these figures that, when the chemical potential exceeds the critical 



3 Different choices of the value of G\ do not change the results in principle, but the scales of parameters 
are changed significantly (see below). 



11 




Figure 2: Condensates a and A as functions of fi for R = 12, Gi — 10, T — (all 
quantities are given in units of A). 




Figure 3: Condensates a and A as functions of fj, for R — 12, Gi = 20, T = (all 
quantities are given in units of A). 

value He, the chiral symmetry is restored and the color symmetry is broken^. 

One may ask, why the diquark condensate A vanishes for large chemical potential /i, 
i.e., why CSC will be destroyed at large densities? (An analogous behavior is seen also 
in Fig. 6 of [23].) In fact, since in the weak coupling QCD the nonvanishing CSC phase 
is present for arbitrary large values of fi, the NJL-results for large values of the chemical 
potential could be an artefact of the model and should be considered with care. 

As is clear from Figs. [H [21 El HI a first order phase transition takes place at the point 
fi = fi c (the first derivative of the themodynamical potential is discontinuous). 

In Fig. [5} the /i — i?-phase portrait of the system at zero temperature is depicted 
for two choices of G\ = 10 and G\ = 20 (1, 2 and 3 denote the symmetric, chiral and 
superconducting phases, respectively). 

For points in the symmetric phase 1, the global minimum of the thermodynamic 

4 It should be mentioned that for R — 3, G\ — 10 our value for // c w 1.75 is close to that from 40J. 
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n 




potential is at a = 0, A = (chiral and color symmetries are unbroken). In the region 
of phase 2, only chiral symmetry is broken and a ^ 0, A = 0. The points in phase 3 
correspond to the formation of the diquark condensate (color superconductivity) and the 
minimum takes place at a = 0, A ^ 0. The dotted line in Fig. [5] denotes the first order 
phase transition, and by the solid line the second order phase transition is depicted. The 
bold point denotes the tricritical point between transitions of the first and second orders. 
As was pointed out in |40j, the vertical line at T = between phases 1 and 2 implies 
that the lowest eigenvalue of ujq is non-vanishing. In the vicinity of this line, a second 
order phase transition takes place, i.e., in this case, the curvature R behaves exactly in 
the same way as the temperature T does in the flat space. Let us mention, that, as in the 
flat space (see, [22]), at T = 0, no mixed phase with both condensates different from zero 
takes place. Comparison of the plots for values G\ = 10 and G\ = 20 demonstrates no 
significant modification of the form of the phase portrait made by the different choices of 
the coupling constant. However, with growing value of G\, the critical values of R and fi 
increase. 

Moreover, the oscillation effect clearly visible in the phase curves of Fig. [5] should be 
mentioned. This may be explained by the discreteness of the fermion energy levels (|2"oT) 
in the static gravitational field. This effect may be compared to the similar effect in the 
magnetic field H, where fermion levels are also discrete (the Landau levels). The corre- 
sponding magnetic oscillations of the \x — y^if -phase portrait in dense cold quark matter 
with four- fermion interactions were found in paper [llj, where it was demonstrated that in 
the massless case, such a phase structure leads unavoidably to the standard van Alphen-de 
Haas magnetic oscillations of some thermodynamical quantities, including magnetization, 
pressure and particle densitj@. In addition, we considered also phase transitions at finite 

5 On the other hand, since these oscillations are related to the restoration of color symmetry for large 
/i, they seem to be here rather an artefact of the NJL model. 
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10 20 30 40 10 20 30 



Figure 5: The phase portrait at T=0 for G\ = 10 (left picture) and G\ = 20 (right 
picture). Dotted (solid) lines denote first (second) order phase transitions. The bold 
point denotes a tricritical point. The numbers 1,2 and 3 designate the symmetric, chiral 
and superconducting phases, respectively. 

temperatures. In Fig. El \i — R- and T — ft- phase portraits are depicted. No mixed phase 
with simultaneous nonzero values of a ^ 0, A ^0 can be seen in these plots. 

It is clear from Figj6] that with growing temperature both the chiral and color sym- 
metries are restored. The similarity of plots in R — \i and \i — T axes leads one to the 
conclusion that the parameters of curvature R and temperature T play similar roles in 
restoring the symmetries of the system. For comparison, we also studied the phase por- 
trait for the case of vanishing curvature, R = 0. The obtained result turns out to be very 
similar to the right picture in Fig. El In this case, with chemical potential growing from 
zero, the critical temperature for the transition between phases 1 and 2 decreases from a 
value slightly higher than 1.5. 




Figure 6: The phase portraits at T=0.35 (left picture) and at R=3 (right picture), G\ = 
10. 
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7 Conclusion 



By using a static model of curved space-time (the Einstein universe), we have derived a 
nonperturbative expression for the effective potential of the Nambu-Jona-Lasinio model 
in the mean field approximation. The corresponding effects of gravitation were exactly 
taken into account. This allowed us to consider phase transitions in the vicinity of the 
critical points. In particular, we have studied the influence of the chemical potential, the 
temperature and the curvature on the chiral and color superconducting phase transitions. 
Moreover, an oscillation effect of the phase curves was found, which may be explained by 
the discreteness of the fermion energy levels in the static gravitational field. 

The analysis made in this paper demonstrates that the space curvature R in certain 
cases plays an analogous role as the temperature T does in the flat space, leading to a 
second order phase transitions and the restoration of chiral symmetry. On the other hand, 
at certain values of the chemical potential (and a sufficiently large coupling constant) a 
superconducting phase is formed, and the color symmetry is broken. However, as might be 
expected, the chiral and color symmetries become restored at high temperatures and high 
curvatures. Clearly, the simple schematic model of gravitation considered in this work 
may represent, at best, only qualitative features of the diquark formation, which possibly 
might occur in the cores of heavy compact objects like neutron stars. Much further work 
has to be done towards a more quantitative and realistic description of the dynamical 
breaking of chiral and color symmetries in curved space-time. This concerns, in particular, 
the important question, whether the predicted CSC in cold heavy compact objects like 
neutron stars might eventually be destroyed by their associated strong gravitational fields. 

Notice that our investigations were primarily concentrated on the study of CSC in the 
4-dimensional Einstein universe. It should, however, be mentioned that the phase curves 
plotted for space dimensions different from 3 demonstrate a considerable similarity with 
the case studied in this work. 
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